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ABSTRACT 


We have used the Spitzcr satellite to monitor the mid-IR evolution of SN 
1987A over a 5 year period spanning the epochs between days ~ 6000 and 8000 
since the explosion. The supernova (SN) has evolved into a supernova rem- 
nant (SNR) and its radiative output is dominated by the interaction of the SN 
blast wave with the pre-existing equatorial ring (ER). The mid-IR spectrum is 
dominated by emission from ~ 180 I< silicate dust, collisionally-heated by the 
hot X-ray emitting gas with a temperature and density of ~ 5 x 10 6 K and 
~3x 10 4 cm“ 3 , respectively. The mass of the radiating dust is ~ 1.2 x 1()~ 6 M 0 
on day 7554, and scales linearly with IR flux. Comparison of the IR, data with 
the soft X-ray flux derived from Chandra observations shows that the IR-to-Xray 
flux ratio, IRX , is roughly constant with a value of 2.5. Gas-grain collisions 


1 Observational Cosmology Lab., Code 665; NASA Goddard Space Flight Center, Greenbelt, MD 20771, 
U.S.A., e-mail: eli.dwek@nasa.gov 

2 CRESST/UMBC, Code 665, NASA Goddard Space Flight Center, Greenbelt MD, 20771, U.S.A. 

•'* DSM/D AP NIA/ Service d’Astrophysique, CEA/Saelay, F-91191 Gif-sur- Yvette? Patrice.Bouchet@cea.fr 

4 Department of Astronomy and Astrophysics, Pennsylvania State University, 525 Davey Laboratory, 
University Park, PA 16802, U.S.A. 

5 Harvard-Smithsonian, CfA, 60 Garden St., MS-19, Cambridge, MA 02138, U.S.A. 

6 Osservatorio Astronomico di Trieste, Via Tiepolo, 11, Trieste. Italy 

7 Gemini Observatory, Southern Operations Center, c/o AURA, Casilla 603, La Serena, Chile 

^Department of Astronomy, University of Minnesota, 116 Church St., SE, Minneapolis, MN 55455, U.S.A. 

University of Wisconsin, Eaii Claire. WI 54702, U.S.A. 



- 2 - 

therefore dominate the cooling of the shocked gas. The constancy of IRX is 
most consistent with the scenario that very little grain processing or gas cooling 
have occurred throughout this epoch. The shape of the. dust spectrum remained 
unchanged during the observations while the total flux increased by a factor of 
~ 5 with a time dependence of p 0 ' 87±0 - 20 . f being the time since the first en- 
counter between the blast wave and the ER. These observations are consistent 
with the transitioning of the blast wave from free expansion to a Sedov phase as 
it propagates into the main body of the ER, as also suggested by X-ray observa- 
tions. The constant spectral shape of the IR emission provides strong constraints 
on the density and temperature of the shocked gas in which the interaction takes 
place. The IR, spectra also suggest the presence of a secondary population of 
very small, hot (T > 350 K), featureless dust. If these grains spatially coexists 
with the silicates, then they must have shorter lifetimes. The data show slightly 
different rates of increase of their respective fluxes, lending some support to this 
hypothesis. However, the origin of this emission component and the exact nature 
of its relation to the silicate emission is still a major unsolved puzzle. 

Subject headings: ISM: supernova remnants - ISM: individual (SNR 1987A) - 
ISM: interstellar dust Infrared: general X-rays: general 


1. INTRODUCTION 

About 10 years after its explosion on February 23, 1987, supernova (SN) 1987A has 
evolved from a supernova, when its radiative output was dominated by the release of ra- 
dioactive decay energy in the ejecta, into a supernova remnant (SNR), when its radiative 
output became dominated by the interaction of its blast wave with the inner equatorial ring 
(ER). The ER is located at a distance of about 0.7 lyr from the center of the explosion, and 
could have been produced by mass loss from a single rotating supergiant (Heger & Langer 
1998) or by a merger event in a binary system which also formed the two outer rings (Morris 
& Podsiadlowski 2009). The transition from SN to SNR was observed at wavelengths ranging 
from radio to X-rays, and is summarized in Figure 18 in Bouehet et al. (2006). 

In this paper we report on the continuing evolution of the ~ 5 - 30 pm spectrum and 
the 3.6. 4.5, 5.8, 8.0. and 24 pm photometric fluxes from SN 1987A, spanning the ~ 5 year 
period from day ~ 6000 until day ~ 8000 after the explosion. Initial reports and analysis of 
the IR evolution were presented by Bouehet et al. (2004, 2006) and Dwek et al. (2008). In 
§2 we present the IR data obtained by the Spitzer satellite. The evolution of the IR emission 
and the dust composition are presented in §3. In §4 we derive the plasma conditions from 



the IR. observations. The comparison of the IR emission with the X-ray emission, and the 
evolution of their flux ratio are discussed in §5. A brief summary of the paper is presented 
in §6. 


2. OBSERVATIONS AND DATA REDUCTION 


Spitzer’s observations using IRAC, IRS, and MIPS were conducted annually during 
the first two years (2004 2005), and then roughly every 6 months for the remainder of its 
cryogenic mission. SN 1987A was also observed incidentally with IRAC and MIPS by the 
SAGE survey (Meixner et al. 2006) and other projects during 2005 . 

Figures 1 and 2 illustrate the IR AC data at the final observations (~ day 8000). The 
IRAC data are all of good quality, i.e. the data are not adversely affected by any artifacts 
induced by bright sources in the field or in prior observations. Initial tests of generating 
mosaics from the basic calibrated data (BCD) showed no significant difference from the post- 
BCD mosaics. Therefore, aperture photometry was performed on the post-BCD mosaics in 
order to obtain the broad band evolution of SN 1987A at 3.6, 4.5, 5.8 and 8 pa n. These 
measurements and their statistical uncertainties are shown in Figure 3. 

Figures 1 and 2 also illustrate the MIPS 24 and 70 /mi images of SN 1987A at Day 7983. 
There is no evident point source at 70 /mi. As for IRAC, there seemed to be no advantage to 
reprocessing the BCD, thus aperture photometry was performed on the post-BCD mosaics. 
The 24 /mi flux densities are shown in Figure 3. At Day 7983, a 70 pm upper limit, was 
determined by adding an artificial point source to the MIPS image at the location of the SN. 
We started with a source flux that, was clearly visible above the background confusion, and 
then reduced the flux until the source was no longer detected using the FIND routine of the 
IDLASTRO library. The results indicate the <S'(70/mi) < 0.09 Jy. 

A significant complication with the IRAC and MIPS photometry is that the SN is 
not resolved from the nearby Star 2 and Star 3. If the stars arc assumed to have mid-IR 
magnitudes that are equal to the K magnitudes (15.06, 15.80 respectively from Walborn et 
al 1993), then their combined IR flux densities should be 0.41, 0.26, 0.16, 0.09, and 0.01 
m.Jy at 3.6, 4.5 5.8. 8, and 24 /mi, respectively. These assumptions may be accurate for Star 
2, which is identified as a B2 III star (Walborn et al. 1993; Scuderi et al. 1996), however 
Walborn et al. (1993) report both, an IR excess and a variability by at least Am = 1 mag 
for Star 3. The fluxes in Figure 3 are plotted after subtraction of the emission of Stars 2 and 
3 under these assumptions. 


IRS observations were reduced starting with the BCD, and using the basic steps out- 



lined in the IRS Data Handbook. For the low resolution (R ~ 60 — 120) data, we con- 
structed a 2-D background frame using the observations of a specific off source target po- 
sition. IRSCLEANJV1ASK was used to identify and clean rouge pixels in the background 
subtracted BCD frames. The cleaned BCD data were averaged for each nod (positions of the 
source at 1/3 or 2/3 the length of the slit), and spectra were extracted from these super BCDs 
using SPICE. After extraction, the spectra orders were merged to produce the full spectra 
shown in Figure 3 (top right panel). The high resolution (R ~ 600) data were processed 
in a similar fashion. The major difference was that instead of using a background from a 
dedicated off-source target, the background is derived from the average of two pointings that 
flank the SN on opposite sides (15" and 25" away for the SH and LH spectra respectively). 
The position angle of these pointings varied as dictated by the scheduling of the observations, 
but the distance of the pointings from the SN remained constant. For both the high and low 
resolution data, SMART was used to fit and extract line fluxes. 


3. THE EVOLUTION OF THE INFRARED EMISSION AND THE DUST 

COMPOSITION IN THE ER 

Figure 3 depicts the evolution of the IRAC and MIPS broad band photometry (top left 
panel), and the low-resolution IRS spectrum of SNR 1987A (top right panel), as a function 
of time. The relatively strong background may affect the continuum at A > 20 /mi and the 
[Ne III] 15.6 pm and [S III] 18.7, 33.5 pin emission lines. Line emission from SN 1987A 
and its surrounding medium will be presented elsewhere (Arendt et al; in preparation). 
The figures show that the IR flux is rising smoothly in all bands, increasing by a factor of 
~ 3 between days 6190 and 7980. The bottom left panel of the figure shows the spectra 
normalized by least-square fit to the one obtained on day 7554. The spectrum maintained 
a nearly constant shape throughout, the 1800 days of observations. This constancy enables 
the use of any IRAC or 24 pm MIPS band as a proxy of the total integrated IR flux from 
the ER. The bottom right panel in the figure shows the ratio between the integrated IR flux 
and the 24 pm flux from the ER. The ratio is constant with a value of 1.64 to within < 10%. 
The constancy in the spectral shape also places strong constraints on the parameters of the 
hot X-ray gas in which the grains are embedded. 

Figure 4 (left panel) shows that the mean spectrum of the ER. can be well fitted by 
astronomical silicate grains radiating at a single temperature of ~ 180 K. The relative ratio 
between the two peaks of the silicate features provide strong constraints on the temperature. 
The mass of the silicate dust is 1.2 x HP 6 Mp on day 7554. and should scale linearly with 
the IR flux for other epochs (see Table 1). The figure also shows the presence of an IR 



excess at wavelengths between ~ 5 — 8 /mi which cannot be attributed to emission from the 
silicate grains. This excess emission was first detected by Bouchet et al. (2006), but was not 
further pursued because of the shorter integration time and larger noise of that early Spitzes 
spectrum. An unidentified K-band (2.2 pm) continuum at a level of ~ 0.4 mJy, detected by 
Kjasr et al. (2007), may be part of this emission component. This emission component, now 
clearly present in all epochs, can be approximated in the 5 to 8 pm wavelength region by 
a u a power law with a spectral index of a w —1.8. This rules out thermal bremsstrahlung, 
which has a well defined spectral index of a ~ —0.1 in this frequency regime, as a possible 
source. Synchrotron radiation can also be ruled as a source of this emission component. The 
observed radio synchrotron emission (Zanardo et al. 2010) has a spectral index of ss —0.8, 
mid its extrapolation from a value of 300 mJv at 1.4 GHz on day 7000 will make a negligible 
contribution to the 8 pm emission around that epoch. We therefore conclude that this 
emission arises from an additional hot dust component that is tightly correlated with the 
silicate emission. Its spectral shape should be featureless at longer wavelengths as to not 
“spoil” the fit to the observed 9.7 and 18 pm silicate emission features. The right panel of 
the figure shows the fit of a blackbody spectrum, modified by a A -2 emissivity law and a 
temperature of ~ 350 K, to this component. 

The modified blackbody is of course an idealized representation of this secondary dust 
component. Figure 5 shows a two-component fit to the IR spectrum when the secondary 
component was taken to consist of either carbon, pure iron, magnetite (Fe 3 0 4 ), or iron sulfide 
(FeS) grains. All components, combined with the main silicate dust component, provided 
acceptable fits to the overall IR spectrum of the ER. The temperature of the secondary dust 
component is generally significantly higher than that of the silicate grains, ranging from 
~ 370 K for pure iron grains to ~ 460 K for carbon grains. Some other possible secondary 
dust components, in particular silicates and aluminum oxides, were ruled out because of the 
presence of strong emission features at wavelengths above ~ 8 pm. 

The problem of accounting for the widely different temperatures of the two dust com- 
ponents, even though they may be collision ally- heated by the same gas, will be addressed in 
FI below. 


4. PLASMA CONDITIONS AS DERIVED FROM IR OBSERVATIONS 

Figure 6 shows the plasma densities required to collisionally heat the different grain 

types to their respective temperatures given in Figure 5 as a function of grain radius. The 
grains are assumed to be embedded in an ionized plasma with a temperature of 5 x 10 4 * 6 K, 
reflecting the average temperature of the soft X-ray component (see §5 below). If we require 



both dust components to arise from the same gas, then only the carbon grains support a 
range of plasma densities (enclosed between the two horizontal dashed lines) that overlap 
with that of the silicate grains. So the IR spectrum of the ER. is consistent with the emission 
from a population of small ( a < 0.05 pm) carbon grains intermixed with a population 
of large (a > 0.20 pm) silicate grains embedded in a plasma with a density n e as n. H « 
(2 - 4) x 10 4 cm ' 3 . Requiring a common environment for the two dust components rules out 
the other dust compositions (Fe, Fe 3 () 4 , and FeS) as viable dust candidates. Alternatively, 
the emission could arise from these grains, provided they reside in a significantly denser 
phase (n e > 2 x 10 5 cm -3 ) of the X-ray emitting gas. This phase then has to be tightly 
correlated with the lower-density one in order to explain the coupling between the two dust 
emission components . 

The potential presence of carbon with the silicate dust in the ER. is surprising, defying 
the common paradigm that either carbonaceous or silicate dust can form in the wind of 
the progenitor, depending on the C/O abundance ratio in the outflow. The observations 
therefore suggest either that CO formation did not exhaust all the carbon in the outflow, or 
that the co-existence of the two dust components is a manifestation of a binary origin of the 
ER. The progenitor of SN 1987 A could have shared a common envelope with a less massive 
carbon-rich star. The interaction between them could have led to the formation of the ER. 
as well as the two outer rings observed by the Hubble to be in the vicinity of the SN (Morris 
& Podsiadlowski 2009). 


5. THE EVOLUTION OF THE INFRARED and X-ray EMISSION AND 

THEIR RATIO 


A rise in the IR emission is expected as the shock sweeps up more dust during its 
propagation into the dusty circuinstellar medium. The rate of increase of the IR emission 
depends on the morphology of the medium into which the shock is expanding, the expansion 
rate of the shock, and on the possibility that dust may be destroyed in the postshock gas. 
The expected increase of the IR emission with time is schematically illustrated in Figure 7 
for the case in which the grain lifetime is sufficiently long that it. can be ignored ( r sput > t', 
where t' is the residence time of the dust in the hot gas; left column), and for the case in 
which the grain lifetime is sufficiently short that an equilibrium is established between the 
rate of grain destruction and the rate at which the dust is replenished by the expanding 
shock ( r sput < t'\ right column). For both eases the figure depicts a shock expanding into 
a 1 -dimensional protrusion (top), a 2-D disk, or ring (middle), and a homogeneous sphere 
(bottom). The dependence of the volume of the swept-up gas is presented for two different 



types of blast waves: (1) a freely expanding one, with a constant velocity v s ; and (2) a blast 
wave that entered the Sedov phase, for which v a ~ if~ 3 f 5 . 

When r sput > if the intensity of the IR emission scales as the mass of swept up dust 
which, for a constant dust-to-gas mass ratio, simply scales as the volume of the shocked gas. 
For a freely-expanding blast wave, the volume of the shocked gas, 1//, will increase as t', 
t 72 , and F 3 , for the three morphological scenarios, respectively. During the Sedov phase the 
volume of the shocked gas, V s , will be increasing of increasing more slowly as d 2 / 5 , d 4</5 , and 
d 6 / 5 , respectively. 

When T sput < d the dust column density through the shock will eventually reach a 
constant value, roughly equal to v s x r sput , where r sput is the grain destruction time in the 
postshock gas. The IR intensity behind a freely-expanding shock will therefore reach a 
constant value if the shock is expanding into a one dimensional protrusion, and will increase 
as d or d 2 when it expands into a ring or sphere, respectively. During the Sedov phase the 
volume of the shocked gas, V 8 , will actually be decreasing as d -3/5 , d -1 / 5 , or increasing more 
slowly as d 1/5 , for the three morphological scenarios, respectively. 

Figure 8 depicts the evolution of the IRAC and 24 /mi MIPS fluxes. The fluxes were 
fitted by a power-law fit. of the form ( t — t 0 ) a , where t is the time since the explosion, and a 
and t 0 are free parameters. The value of t 0 can be interpreted as the time when the shock 
first, encountered the ring. The values of {a,t 0 } were found to be {0.6, 5700 d} and {0.87, 
5500 d} for the IRAC and MIPS fluxes, respectively. The relatively large value of a suggests 
that grain destruction is unimportant. More specifically, the 24 n m flux increases at a rate 
that is consistent with a Sedov blast wave expanding into a circumstellar disk/ring. 

The value of t Q is consistent with the time inferred from X-ray observations that showed 
a deceleration of the radial expansion of the SN around day 6000 (Racusin et al. 2009), and 
the upturn in the soft X-ray luminosity that occurred around day 6200 (Park et al. 2006). 
The IR and X-ray observations therefore offer a consistent scenario in which both fluxes 
originate from the dust and gas shocked by the blast wave that has started to penetrate the 
main body of the ER around day 6000 after the explosion. 


5.1. The Evolution of IRX 

Analysis of the Chandra X-ray observations of the ER show that the spectrum is 
generated by two distinct emission components: a soft component with temperatures ~ 
(0.3— 0.6) keV, and a harder component characterized by temperatures between ~ (2—5) keV 
(Zhekov et al. 2009: Park et al. 2006). The emission from the soft component arises from the 



gas heated by the shock that is transmitted into the ER, whereas the hard component arises 
from the lower density gas interior to the ER that is shocked once by the advancing SN blast 
wave, and a second time by the shock reflected from the ER. We assume here that the IR. 
emission arises from the soft X-ray component generated by the shock that is penetrating 
the ER. 

Figure 9 (left panel) depicts the evolution of the soft X-rays and IR fluxes, and their 
ratio, IRX , as a function of time. The IRS observations only cover the ~ 6804 — 8000 d 
epoch, so we used the MIPS 24 pm observations that cover the additional ~ 6190 — 6800 d 
epoch as a proxy for the integrated 5-30 pm emission. The figure also depicts the day 6067 
datum point, for which we used the scaled T-ReCS 10.4 pm observations (Bouchet et, al. 
2006) as a proxy for the integrated IR flux. 

The evolution of the radius of the soft X-ray emission component shows a distinct flat- 
tening in the slope around day ~ 6000, which was interpreted as evidence for the transition 
of the blast wave from the free-expansion to the Sedov phase as it enters the main body 
of the ER (Park et al. 2005; Racusin et al. 2009). The T-ReCS datum suggests that this 
transition may have also manifested itself in the IR, emission. 

The evolution of IRX is shown in more detail in the right panel of the figure. It has a 
roughly constant value of ~ 2.5 throughout the entire ~ 6000 — 8000 d epoch. An exception 
is the 6190 datum point. We believe that the anomalously large value of IRX at this point 
is caused by a relative dip in the X-ray flux starting near day 5800 and ending with the rise 
in the soft X-rays near day 6200. 

A constant value of IRX requires the IR and X-ray fluxes to evolve at the same rate. 
This requires the dust-to-gas mass ratio to be uniform throughout the medium into which 
the shock is expanding, and that neither grain destruction nor gas cooling affect the IR and 
X-ray emission from the shocked gas. Alternatively, the constancy in the ratio could be the 
result that both processes, grain destruction and gas cooling, are important, and that the 
IR and X-ray emission have attained a steady state in which the losses to grain destruction 
and cooling are balanced by the influx of new material. The rate of increase in the IR, fluxes 
suggest that the second scenario is a-priori not likely to be the cause for the constancy in 
IRX. In flit' following we examine the physical condition of the shocked gas that can lead 
to a constant IRX in more detail. 



9 - 


5.2. Grain Destruction in the Hot Gas 

The lifetime of a dust grain of radius a moving at a high velocity through a hot gas with 
temperatures between ~ 1() 6 — ID 7 K is given by (Nozawa et al. 2006; Dwek et al. 2008): 

/ \ . a(am ) , , 

Tsputiyr) = A sput Mcm — 3) (1) 

where A sput = 3.3 x 10°, 9.0 x 10 5 , 5.6 x 10 5 , 2.8 x 10 5 ,and 3.9 x 10°, respectively, for silicate, 
carbon, iron, FeS, and Fe 3 0 4 dust, calculated for grains moving at velocities > 500 km s _1 
through a 5 x 10 6 K gas. Figure 10 depicts the sputtering lifetime of the dust as a function of 
grain radii and gas density. The figure shows that the lifetime of the silicate grains is between 
~ 4 and ~ 15 yr, depending on the combination of grain size and gas density, whereas that 
of the carbon dust is between ~ 0.4 and ~ 1 yr. 


5.3. Cooling of a Shock-Heated Dusty Plasma 

Figure 11 depicts the equilibrium atomic cooling rate, A (T e ), of a gas with ER compo- 
sition (thick line, left panel) as a function of electron temperature, T e . The ER. abundances 
were taken from Mattila et al. (2010), and were derived from AAT and VLT observations 
of the supernova between days ~ 1400 and ~ 5000 after the explosion. The thin lines in 
the figure show the cooling of the plasma via gas-grain collisions, assuming a pure silicate 
dust composition. The silicate dust-to-gas mass ratio was assumed to be 0.0029, half the 
solar value. The different curves assume that all the dust is represented by grains with the 
indicated radii. Dust cooling dominates the cooling via atomic transitions at temperatures 
above ~ 10 6 K, depending on grain size. The value of IRX is the ratio of these two cooling 
rates. From the fact that its initial value is only ~ 2.5 we can conclude that the grain sizes 
in the ER must be relatively large, with typical radii of > 0.2 /on, consistent with the sizes 
required to heat the dust to a temperature of ~ 180 K (see Fig. 6). Using this lower limit 
on the grain radii in eq. (1), gives a postshock density of about (2 — 4) x 10 4 cm -3 , where 
we assumed an upper limit of 0.5 /mi on the grain radii. 

The equilibrium cooling may not represent the actual cooling of a shocked gas, which 
may be far from ionization equilibrium. The right panel in Figure 11 depicts the temperature 
profile behind shocks of different velocities expanding into a homogeneous dusty gas with 
a preshock hydrogen density of 10 4 cm -3 and ER composition. The shock models include 
the effects on non-equilibrium ionization (NEI), and the effect of dust cooling. The ~ 
600 — 700 km s^ 1 shocks well represent the 3 — 7 keV temperature of the soft X-ray emitting 
plasma (Park et al. 2006; Zhekov et al. 2009). The figure shows that the cooling time of 
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this shocked gas is about 12 — 20 yr for a presliock gas density of n H = 10 4 cm -3 . The IR. 
observations suggest a postshoek density of ~ (2 — 4) x 10 4 cur 3 , giving a gas cooling time 
of ~ 20 — 40 yr. 


5.4. Scenarios for the Constancy of IRX 

The results presented above show that the timescales for grain destruction and gas 
cooling are sufficiently long that neither processes has significantly affected the IR. or X-ray 
emissions, or their flux ratio, IRX. The sputtering timescale suggests that grain destruction 
may become important only at day ~ 6000 + 3200 ~ 9200 after the explosion, and the 
cooling time suggests that the X-ray emission may not be affected until t ss 30 yr. 

The alternative scenario suggests that the approximate constancy in IRX is the result 
of both processes, grain destruction and gas cooling, having attained a steady-state, in which 
the amount of preshocked gas and dust entering the shock is balancing the amount of dust 
that is destroyed and the gas that has cooled. Such scenario requires the grain destruction 
timescale and gas cooling time to be significantly shorter than ~ 100 d. requiring postshoek 
densities that are more than an order of magnitude larger than those inferred from the 
temperature of the collisionallv- heated silicate dust. 

The original scenario in which the evolution of the IR and X-ray emission from the 
postshoek gas has not yet been affected by the two processes offers therefore a more rea- 
sonable explanation for the constancy of IRX. It relies on the existence of a population of 
relatively large, a > 0.2 pm, silicate grains, embedded in a gas with postshoek densities of 
about (2 — 4) x 10 4 cm -3 . These numbers should be regarded as suggestive values for the 
model parameters. 

However, this simple explanation is complicated by the presence of the secondary dust 
component consisting of much smaller dust grains. Figure 10 shows that the density-radius 
combinations required to heat the carbon dust to ~ 460 K will also destroy the grains in less 
than half a year, and the various iron bearing grains would last for less than ~ 10 d, which 
definitely rules them out as viable candidates for the secondary dust component. 

The short sputtering time for the carbon dust suggests that their column density through 
the shock reached a steady-state between injection and destruction shortly after the shock 
penetrated the ER. Consequently, the IR emission from these dust particles should have 
increased at a different rate than the emission from the silicate grains. This is in contrast to 
the good correlation between the two emission components. 
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Figure 8 shows that the 5 — 8 /xm emission does increase slightly less rapidly than 
the 24 /mi emission. However, the difference in slope is only 1/4 as large as expected if 
the r sput (hot dust) < t' < r sput (silicate). It is therefore possible the ~ 5 - 8 /an emission 
could arise from a yet unknown dust component consisting of grains with either radii or IR 
emissivities (or both) sufficiently small to reach the observed high temperature, and with 
binding energies sufficiently large so that their sputtering lifetime will be ~ 5 yr in spite of 
their small radius. Furthermore, they should have a featureless spectrum in the ~ 4- 40 /an 
wavelength region, in order not to ” spoil" the good fit to the ~ 9 — 20 /on silicate spectrum. 
They could also be embedded in a different, lower density, environment than the silicate 
dust, which will increase their lifetime. The nature of this secondary dust component is still 
an unresolved aspect of the observations. 


6. SUMMARY 

In this paper we presented the mid-IR evolution of SN 1987 A over a 5 year period 
spanning the epochs between days ~ 6000 and 8000 since the explosion. Its radiative output 
during this epoch is dominated by the interaction of the SN blast wave with the pre-existing 
equatorial ring (ER). The main results of this paper can be briefly summarized as follows: 


1. The ~ 8 - 30 /xm mid-IR spectrum is dominated by emission from ~ 180 Iv silicate 
dust, collisionally-heated by the hot X-ray emitting gas with a temperature and density 
of ~ 5 x 10 6 K and ~ (2 — 4) x 10 4 chi'" 3 4 , respectively. The mass of the radiating dust 
is ~ 1.2 x 10~" 6 Mq on day 7554, and scales linearly with IR flux. 

2. A secondary emission component dominates the spectrum in the ~ 5 - 8 /xm region. 
Its intensity and spectral shape rule out any possible gas or synchrotron emission 
mechanism as the source of this emission. It must therefore attributed to a secondary 
dust component radiating at temperatures above ~ 350 K. 

3. The overall shape of the ~ 5 - 40 /xm dust spectrum remained unchanged during the 
observations, suggesting that the density and temperature of the soft X-ray emitting 
gas have not significantly changed during the more than 5 years of IR observations. 
The constancy in the spectral shape of the IR. emission also suggests that the mass 
ratio of the silicate to the secondary dust component remained roughly constant during 
this epoch 

4. The evolution of the IRAC and MIPS fluxes can be described by a power law in time 
since the first shock-ER encounter. The silicate emission increases as R 8 ' , consistent X- 
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ray observations suggesting that the blast wave has transitioned from a free expansion 
to the Sedov phase, and is now expanding into the main body of the ER. 

5. The infrared-to-X-ray flux ratio, IRX, is constant with a value of ~ 2.5 throughout this 
epoch. The magnitude of IRX shows that the cooling of the shocked gas is dominated 
by 111 emission from the collisionally-lieated dust with radii > 0.2 fan, and that a 
significant fraction of the refractory elements in the ER should be depleted onto dust. 

6. The constancy of IRX is consistent with the premise that neither grain destruction by 
sputtering nor cooling of the shocked gas have played a significant role in the evolution 
of the 1R and X-ray emission. 

7. The presence of a secondary dust component, radiating at significantly higher temper- 
atures than the silicate dust, suggests that the grain radii or IR emissivities of this 
component must be significantly smaller than those of the silicates. Their sputtering 
lifetime could therefore be significantly shorter than that of the silicate grains and their 
evolution distinctly different from that of the silicates, especially around the ~ 6100 d 
time interval. The constant mass ratio between these two components during this 
epoch cannot be explained by the common dust candidates considered here for the 
secondary dust component. 


We benefitted from useful comments by Dick McCray and Svetozar Zhekov. This work 
is based in part on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology under a contract with 
NASA. Support for this work was provided by NASA. 



13 

REFERENCES 

Bouchet, P.. De Buizer, .J. M., Suntzeff, N. B., et al. 2004, ApJ, 611, 394 

Bouchet, P., Dwek, E,, Danziger, J., et al. 2006, ApJ, 650, 212 

Dwek, E., Arendt, R. G., Bouchet, P., et al. 2008, ApJ, 676, 1029 

Heger, A. & Langer, N. 1998, A&A, 334, 210 

Kjaer, K., Leibundgut, B., Fransson, C., et al. 2007, A&A, 471, 617 

Mattila, S., Lundqvist, P., Groeningsson, P., et al. 2010, ArXiv e-prints 

Meixner, M., Gordon, K. D., Indebetouw, R... et al. 2006, AJ, 132, 2268 

Morris, T. & Podsiadlowski, P. 2009, MNRAS, 399, 515 

Nozawa, T., Kozasa, T., & Habe, A. 2006, ApJ, 648, 435 

Park, S„ Zhekov, S. A., Burrows, D. N., et al. 2006, ApJ, 646, 1001 

Park, S., Zhekov, S. A., Burrows, D. N„ k McCray, R. 2005, ApJ, 634, L73 

Racusin, J. L., Park, S., Zhekov, S., et al. 2009. ApJ, 703, 1752 

Scuderi, S., Panagia, N., Gilmozzi, R., Challis, P. M., k Kirshner, R. P. 1996, ApJ, 465, 956 
Walborn, N. R., Phillips, M. M„ Walker, A. R., & Elias, J. H. 1993, PASP, 105, 1240 
Zauardo. G., St avelev- Smith, L., Ball, L., et al. 2010. ApJ, 710, 1515 
Zhekov, S. A., McCray, R„ Dewey, D., et al. 2009, ApJ, 692, 1190 


This preprint was prepared with the AAS RTjrX macros v5.2. 















W) ( J y)' vl v (erg c] 


- 16 




6000 6500 7000 7500 8000 

time (day) 


10 

X (iura) 




X (/xm) time (day) 


Fig. 3. Top row: The evolution of the IRAC and MIPS 24 pm luminosities (left panel) and the specific 
intensity of the low-resolution spectra (right panel) of the equatorial ring. Bottom row: Spectra normalized 
to the one obtained on day 7554 (left panel). Because of the constancy in the shape of the spectrum, any of 
the III AC or MIPS photometric fluxes can be used as a proxy for the integrated IP spectrum from the ER. 
The right panel depicts the almost constant ratio between the IR.S integrated over the 5-30 pm spectrum to 
th(' 24 pm luminosity. 
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Fig. 4. Left panel: Single composition fit, to the IR. spectrum of the ER. The fit corresponds to the 
spectrum of astronomical silicate grains radiating at 180 K. The excess 5-8 gm emission evolves at the 
same rate as the silicate emission (see Fig. 3) strongly suggesting that it is generated by an additional 
dust component in the ER. Right panel: The addition of any dust component with a smooth spectrum 
which will not “spoil” the fit to the 9.7 and 18 pm silicate features is a viable candidate for this second dust 
component. Here the spectrum is fit with a modified blackbody with a A -2 emissivity law and a temperature 
of 346 K. 
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Fig. 5. Two-composition fits to the observed infrared spectrum of the ER. indicated by the violet curve. 
An additional hot component consisting of either carbon, metallic iron, iron oxide, or iron sulfite grains 
(red curves) provides a good fit to the 5 — 8 pm segment of the spectrum. The blue curve represents the 
spectrum of the silicate grains, and the green curve represents the sum of the emission from the silicate and 
the secondary dust component. 
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Fig. 6. The density required to collisionally heat the different grain species to their respective tempera- 
tures given in Figure 5 as a function of grain radius. Details in text. 
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Fig. 7.~ A schematic view of the evolution of the volume of shocked swept-up dust for different geometries 
of the cireumstellar medium: a linear protrusion (top row); a hollow disk (middle row); a uniform sphere 
(bottom row). The amount of radiating hot dust is depicted as a dark-shaded region behind the advancing 
shock. The left column depicts the evolution when grain destruction is insignificant, i.e. when T sput > t, 
where t is the time when the shock first encounters the cireumstellar medium. The mass of shock-heated 
dust increases in proportion to the swept up gas. The right column depicts the volume of the shock-heated 
dust when r;,p„, > t. The amount of radiating dust is then constant as the increase in the mass of the 
swept-up dust is balanced by its destruction. Each scenario exhibits a different dependence on time, which 
also depends on tin 1 expansion law of the blast, wave: Vjr is the volume swept up for a freely-expanding 
blast wave (v e ~ t°), and V„ is the volume swept if the blast wave entered the Sedov phase of its evolution 
(v„ ~ t~V 5 ). 
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Fig. 8. The evolution of the (renormalized) IRAC and MIPS 24 yum intensities as a function of t' = t-to, 

where to is the time at which the SN blast wave first encountered the ER. The value of to is equal to 5700 d 
for the IRAC fluxes, and equal to 5500 d for the 24 yum flux. The total IR flux emitted by the silicate dust 
(represented by the evolution of the 24 yum luminosity) exhibits a steeper rise than that of the secondary 
dust component, represented by the IRAC luminosities. The rise in the IR spectrum is consistent with a 
Sedov blast wave expanding into a disk-like medium, with negligible grain destruction see Figure 7. 
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Fig. 9. Left panel: The evolution of the IR flux, represented by the 24 yum MIPS (filled blue diamonds) 
and 10.4 yum T-ReCS (open blue diamond) observations, the soft X-ray flux (derived from the Chandra 
observations), and IRX, the IR-to-X-ray flux ratio. Infrared and X-ray fluxes are given in units of 10~ 12 erg 
s _1 cm -2 . IRX is dimensionless. Right panel: Details of the evolution of IRX. Solid blue diamonds and 
the open green diamond are derived using the MIPS 24 /urn and the T-ReCS 10.4 yum fluxes as proxies for 
tire integrated IR spectrum. The IRX derived using the integrated IR fluxes is shown as filled red diamonds. 
The data are consistent with a constant value of IRX over the entire ~ 6000 — 8000 d epoch. 
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Fig. 10. - The sputtering lifetime of the different grain species calculated for the same gas density- grain 

radii combinations required to heat the different dust species to their observed temperature. Dust species 
are color coded as in Figure 6. Within the permitted range of silicate grains, the straight lines indicate the 
loci of constant sputtering lifetimes T sput — 3, 6, 9 (bold), 12 yr. For carbon dust the lines indicate r sput = 
0.4. 0.6, 0.8 (bold), 1.0 yr. For iron oxide dust the lines indicate r spMt = 0.02, 0.03 (bold), 0.04, 0.05 yr. 




Fig. 11. — Left panel: The equilibrium cooling curve for a plasma with ER. abundances (thick line). 
The thin solid lines represent the gas cooling via dust-gas collisions. All curves assume that all the dust 
is represented by silicate grains with radii indicated in the figure. Right panel: The evolution of the 
temperature behind shocks of different velocities, expanding through a homogeneous dusty medium with ER 
abundances and a preshock hydrogen density of 1 x 10 4 cm' 3 . 
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Table 1. Observed X-ray and Infrared Fluxes From SN 1987 A 1 


day 1 2 

X-ray flux 3 

IR, flux 4 5 

IRX 5 

6067 

(1.21 ±0.18) x 10- 12 

(3.15 ±0.46) x lO'" 12 

2.6 ± 0.5 

6184 

(1.40 ±0.19) x 10“ 12 

(5.38 ± 0.37) x nr 12 

3.8 ± 0.6 

6552 

(2.72 ±1.0) x nr 12 

(7.43 ±0.38) x ltr 12 

2.7 ±1.0 

6734 

(3.47 ±0.9) x 10- 12 

(8.52 ± 0.40) x It)- 12 

2.4 ±0.7 

6829 

(3.62 ± 0.79) x 10“ 12 

(9.07 ± 0.38) x 10“ 12 

2.5 ± 0.6 

7159 

(4.43 ± 0.66) x 10- 12 

(1.12 ±0.04) x 10~ n 

2.5 ± 0.4 

7310 

(4.84 ± 0.90) x 10- 12 

(1.22 ±0.04) x 10- 11 

2.5 ±0.5 

7490 

(5.10 ±0.65) x 10~ 12 

(1.32 ±0.04) x 10~ n 

2.6 ± 0.3 

7690 

(5.83 ± 1.18) x 10“ 12 

(1.43 ±0.04) x 10~ n 

2.5 ± 0.5 

7983 

(5.79 ± 1.00) x 10“ 12 

(1.55 ±0.04) x It)- 11 

2.7 ±0.5 


1 A11 fluxes are in units of erg cm -2 s -1 . 

2 Measured since the explosion. 

3 Soft X-ray flux in the 0.5-2. 0 keV band, interpolated to the 
epochs of the Spitzer / MIPS observations and corrected for an ex- 
tinction column density of Nh = 2.35 x 10 21 cm' 2 . 

4 IR fluxes on days 6184 to 7983 were derived from the MIPS 
24 /iin data using the conversion factor of 1.64 (see Fig. 3). The 
flux on day 6067 was derived using the 10.4 pm T-ReCS flux as a 
proxy for the integrated 5-30 pm spectrum. 

5 The ratio of the IR to soft X-ray flux from the ER. 



